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Abstract 

Several investigators of the Apollo lunar experiments have 
observed gradual increases in the mean temperatures recorded by 
various surface thermometers. Similar effects have been noticed 
in the temperatures of the thermometers of the Apollo 15 and 17 
Heat Flow Experiments. This report discusses an analysis of the 
long term temperature histories of the heat flow experiment thermo- 
meters. These data show that no change in mean surface temperature 
at the Apollo 15 and 17 sites has occurred, and suggest that the slow 
increase in "mean" temperatures of thermometers in the electronics 
housing are due to changes in radiative properties of the housing's 
surfaces. 


Note- The Technical Officer for this Contract is Mr. Wilbert F. Eichelman 
TE6 Lyndon B. Johnson Space Center, Houston, Texas 77058 
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Introduction 

Accurate sets of thermometers were placed on the moon as 
part of the Lunar Heat Flow Experiments. One experiment was installed 
at the Apollo 15 site at Rima Hadley, in July 1971, and one at the Apollo 
17 site in Taurus Littrow, in December 1972. All of these thermometers 
have been returning data to earth since they were installed on the moon, 
and provide us with records of temperature variation over a 3. 5 and a 2 
year period. Some of the Heat Flow Experiment thermometers are above 
the lunar surface and assume temperatures that satisfy a heat balance be- 
tween impinging thermal radiation from the sun and the lunar surface, and 
that lost by radiation to space. These thermometers experience very large 
variations in temperature throughout a lunation and during eclipses (Ref- 
erences 1 and 2). The amplitudes of these monthly variations show a 
strong annual modulation and additional weaker modulations over much 
longer times. In this report we describe long term variations that have 
been observed and compare them with the variations expected. 

One of the Heat Flow Experiments' sensors is inside the elec- 
tronics box housing. The temperatures of this sensor show an annual 
variation and a slow increase in mean temperature with time. 
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Measurements ; 

Thermometers used for the lunar heat flow measurements 
are contained in slender probes placed in predrilled holes in the lunar soil 
and in the cables connecting the probes to the electronics unit. Each 
probe contains eight platinum resistance thermometers. There are four 
thermocouple junctions in each cable, which are located at distances of 
approximately 0, 0. 65, 1. 15 and 1. 65 m from the topmost platinum thermo- 
meter in the probe, see Figure 1. The electronics are contained in a 
thermally <;ontrolled housing. Another platinum resistance thermometer, 
the reference thermometer, is attached to the radiator plate of this housing. 
The reference junctions of all of the thermocouples are thermally connected 
to the reference thermometer. Thermocouple no. 1 is inside the topmost 
platinum resistance thermometer in the probe; thus, this platinum thermo- 
meter provides a second reference. Table 1 (from Reference 3) shows 
the absolute accuracy and range of these thermometers. 

Figure 2 shows the emplacement of the probes at the Apollo 
15 and Apollo 17 sites. At Apollo 15 the holes could be drilled to only half 
the intended depth and, consequently, many of the thermocouples in the 
cable were left exposed above the surface. At probe 1 thermocouple no. 4 
is inside the portion of the borestem above the surface. At Apollo 17 the 
holes were drilled to the desired depth and only thermocouple no. 2 is 


exposed above the surface. 



In normal operation, the temperatures of the heat flow experi- 
ment thermometers are transmitted to earth every 7. 25 minutes. Typi- 
cal thermocouple temperature variations during a lunation cycle are 
shown in the inset in Figure 3. For each sensor two temperatures are 
selected from a lunation cycle and used to examine the long term 
variations: 

1) The maximum temperature near lunar noon, and 

2) the temperature at a prescribed time after lunar sunset 

which we call the "presunrise temperature". 

The prescribed times used are 14 days for Apollo 17 and 15 days for 
Apollo 15, which are just before lunar sunrise and, consequently, these 
temperatures are near the minimum values for each lunation. 

Results of Thermocouple Measurements Above the Lunar Surface 

Temperature variations for three typical thermocouples are 
shown in Figure 3. Daytime maxima for these thermocouples show the 
dominating character of the annual component of the insolation. The 
modulating envelopes of the annual variations are, in general, different 
for each thermocouple, and depend on the orientation of the cable that 
surrounds the thermocouple. The effect of cable orientation can be seen 
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in the temperature variations during a lunation (see Figure 3, inset). 

The individual thermocouples reach their maximum temperatures at 
different times in the lunation. The cable element around thermocouple 
no. Z, at A- 17, is oriented nearly north- south and has a nearly symmetric 
curve relative to lunar noon. Thermocouple no. 4 at A- 15 is in the bore- 
stem projecting above the surface. The borestem tilts slightly toward the 
east. This orientation results in a relatively flat curve that peaks early 
in the afternoon.’ The fact that the temperatures of thermocouple no. 2, 
at A- 15, peak in the afternoon and are substantially colder than the other 
thermometers in the lunar morning indicates that the cable axis orien- 
tation at this junction is roughly northwest- southeast and is elevated to the 
southeast. 

Presunrise temperatures of all thermocouples are much more 
constant with time although there is a small annual variation of about half 
a degree peak to peak. Mean presunrise temperatures show no detectable 
drift, except for an abrupt increase in Apollo 15 temperatures between 
lunation 18 and 19. This single discontinuity in the Apollo 15 presunrise 
temperature curves probably results from a spurious change in extraneous 
EMF in the thermocouple circtaitry. Otherwise, the data indicate that the 
thermocouples are quite stable at night. 



Factors Affecting Long-Term Temperature Variations of the Lunar Surface 


1) The long-term variations in maximum surface temperature 
of a level surface element on the moon depend on the time variation of 
insolation and the lunar latitude of the element. The Appendix gives a 
derivation of solar flux at a fixed point on a smooth, spherical moon as a 
function of time. Figure 4B shows the lunation maxima of this function 

for the periods from July 1971 to October 1974 at Rima Hadley and from 
January 1973 to October 1974 at Taurus Littrow. The insolation at Taurus 
Littrow is more intense because that site is nearer the equator. The 
strong annual component results from the eccentricity of the earth's orbit. 
The amplitude modulation of this annual component results from the pre- 
cession of the moon's spin axis with a period of 18. 6 years. 

2 ) The small variation in presunrise temperature depends on 
the annual variation in total flux during the lunar day and the thermal inertia 
of the lunar regolith. There are two factors involved in determining total 
flux; The variation of insolation intensity and the variation of the length of 
the insolation period. In Figure 4A, the variation in the period from sunrise 
to sunset, based on Equation 1 1 of the Appendix, is shown and compared 
with the observed period determined from the thermocouples. The fractional 
'variation in the length of the insolation period is much smaller than that of 
the insolation intensity, and, consequently, has a proportionately smaller 
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effect on the presunrise surface temperature variations. In Figure 5, the 
theoretical variation in presunrise surface temperature is shown based on 
the insolation function and a thermal model of the lunar surface (see 
References 1 and 2). 

Notice that for presunrise temperatures the phase of the theore- 
tical curve lags that of the annual insolation variation by about ir/4 radians. 
This phase shift is expected because, for a semi-infinite solid, the variation 
of surface temperature lags a periodic variation in flux at the surface by 
exactly 7r/-i radians. See, for example. Reference 4, page 76. 

Thermocouple Temperatures ; 


Thermocouples, buried inside a portion of cable exposed 

above the lunar surface, attain temperatures that provide radiative balance 

between iiiipinging radiation from the sun and lunar surface and the cable's 

radiation to space and the lunar surface. This balance is expressed by: 

_4 


F fa. 
cm m cir 


a s sin B 
m ^ cs n 


ir <7 « 


AF a s cos Qc 
cm cs n 

cr « 


( 1 ) 


The notation is the same as that used in Reference 2, except for s , the 

n 

normally incident flux which is a function of time; oc, the incident angle of 

solar radiation on the lunar surface; and /3 , the angle between the cable axis 

and impinging solar radiation. The absorptivities are denoted by sniall a's. 

At lunar nights = 0 and the thermocouple temperature is 
n 
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proportional to the lunar surface temperature . 

4 A/ 

During the lunar day T =s cos oc (1-A)/ct«0‘. With 

m n 

this assumption, the thermocouple temperature can be expressed solely as 
a function of solar insolation and the angles oc and , i. e. , 


T =8 cos oc (K 
c n 


_8inA + K j 
1 cos oc 2 


( 2 ) 


where all constants have been lumped into and K^. This relation shows 
that the variation of thermocouple temperature with variation of is a 
function of the angle of orientation of cable. For a vertically oriented 
cable, for instance, oc = ^ and 


4 

T = s cos oc (K, tan oc + K.) (i) 

c n 1 2 

whereas for a horizontal, north- south oriented cable sin /3 = cos oc at 
lunar noon and 

= 8 cos oc (K, + K,) (4) 

c n 1 2 

Discussion of Measurements 

Pr^yunn£^T«n£^r^;^jr£8_ 

The amplitudes and phase lags of the presunrise temperatures 
show no significant disagreement with those derived from the thermal model 
of the lunar soil (see Figure 5B).The mean values of these temperatures, 
for the exposed thermocouples at both sites, agree well with expected values 
based on measurements of the thermal and radiative properties of the cables 
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and the lunar surface. 

Lunar D aytirne Temperatures 

The long term modulations of the annual variation are due to 
the changes in the angles oc and /3 as the moon precesses. Equations (2), 
(3) and (4) show how the cable orientation determines the shape of this 
modulating envelope for some special cases. Thermocouple 2, at Hadley 
Rille, which is approximately horizontal shows the same kind of modulation 
as the surface temperature (see Figure 5 and Equation 4). Thermocouple 
no. 4, at the same site.- is in the projecting bore stem and thus has an 
approximately vertical orientation. At lunar noon thermocouple no. 4 is in 
equilibrium with the borestem. Writing Equation (3) as 

4 

T = K.s sin oc +K,s cos oc 
c In 2 n 

shows that two modulations, 90* out of phase, are involved. The second 
term on the right defines the effect of heating by radiation from the lunar 
surface, which has an increasing envelope over the period of observation 
being discussed, while the direct solar heating of the cable element defined 
by the first term on the right, has a contracting envelope. The resulting 
envelope will be a sum of these terms, the exact shape depending on the 
mean value of oc and the relative sizes ofK^ andK^. Figure 3 B shows that 
there is only a very small modulation of thermocouple no. 4 temperatures, 
thus the two modulations appear to cancel. 
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A similar analysis, however, cannot explain the variation in 
the maximum temperatures of thermocouple no. 2 at Taurus Littrov. The 
lower temperatures of this sensor near oerihelion 1974 (see Figure com- 
pared with those near perihelion 1973, could possibly be accounted for by 
cable orientation. However, the same reasoning would require the tem- 
peratures near aphelion 1974 to be higher than those near aphelion 1973. 

The observed maximum temperatures, do not conform to this requirement. 

It is not certain that this anomalous behavior is significant, however, since 
the instruniental noise in this sensor near noon makes a irate temperature 
determinations difficult. 

Besides the specific examples already discussed, there are 
some general patterns in the maximum temperature variations of the 
thermocouples which can be explained in terms of precession and cable 
orientations. It can be seen in Figure 5A that the phase of the maximum 
thermocouple temperatures lead maximum surface temperatures by various 
amounts up to 1 lunation. This is also true for all the thermocouples not 
shown in the figures. In fact, the maximum surface temperatures, which 
are in phase with the maximum insolation curve, maximize each year about 
a lunation after perihelion and minimize about a lunation after aphelion. 

This phase lag is caused by the northward advance of the sun near peri- 
helion each year and the corresponding southwai-d recession near aphelion 
due to the precessional motion, which, in addition to increasing the 
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amplitude of the annual variation, shifts the extrema of this variation 
to the right. Several years from now, when the amplitude of the annual 
wave is decreasing, the extrema will be shifted toward the left. The 
thermocouples are not as strongly affected by these precessional effects, 
because they are, in general, not north- south oriented so that the change 
in the angle /d (in Equation 2) is not, in general, as large af the change in 
oc. Thus, the direct heating of the thermocouples tends to maximize and 
minimize nearer to perihelion and aphelion respectively, causing -he 
phase lead seen in the maximum thermocouple temperatures. 

It has also been noted that there is a significant lack of symmetry 
in the modulating envelope of the average maximum temperatures of the 
Apollo 15 thermocouples, which implies a drift in the mean temperature 
of all these sensors. Precessional effects cannot account for a drift in 
mean temperature of the lunar surface, or a thermocouple junction; but, 
the observed drift can be removed by correcting for an apparent error, 
caused by the circuitry. The difference in temperature between junction 1 
and the topmost platinum resistance thermometer in probe 1 near noon at 
both sites is shown in Figure 6B and D. The platinum sensor is assumed to 
be stable, while the relative temperature of junction 1 shows a steady upward 
drift through lunation 18, a sudden jump between lunations 18 and 19, and 
subsequently it is stable. The temperatures of thermocouples 2 and 4 at 
Hadley Rille, corrected for this drift are shown in Figure 5A. A similar 
correction for the temperatures of thermocouple 2 at Taurus Littrow has not 
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been made since it does not significantly change the anomalous behavior of 
this sensor. 

Reference Thermometer Temperatures 

In Figure 6A and C the maximum temperatures and the presunrise 
temperatures of the reference thermometer in the electronics box are shovm 
for both experiments. Except for the anomalously high values in the second 
lunation at each experiment, caused by increased heating of the electronics 
box during conductivity experiments, the maximum temperatures show the 
annual variation superimposed on a gradual upward drift in mean temperature. 

At Apollo 15, this drift seems to be leveling after three years, whereas the 
period of observation at Apollo 17 is not sufficiently long to confirm a similar 
effect there. The similarity in the initial drift rates at the two sites suggests 
that this effect is characteristic rather than spurious. 

The possibility that this drift in mean temperature is caused by a 
drift in calibration of the reference thermometers is considered unlikely for 
the following reasons. First, the reference thermometers are platinum 
resistance bridges similar in construction to those in the probe bodies which 
preflight tests showed to be extremely stable (Kleven et al. , 1970). Next, 
thermocouple 1 has one junction thermally connected to the reference thermo- 
meter and the other inside the uppermost platinum resistance element in the probe. 
This provides direct comparison of the probe and reference bridge. The 
difference in temperature between thermocouple 1 and the uppermost platinum 
resistance thermometer at the Apollo 17 site (see Figure 6B) near noon shows 
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no long-term drift and indicates that the calibration of the reference thermo- 
meter in the Apollo 17 experiment is stable. Therefore, the upward drift 
shown in Figure 6A is a real change in mean temperature of the electronics 
box. At Hadley Rille (Apollo 15), the same stability is evident after lunation 19 
(see Figure 6D), indicating that the slow increase in mean reference thermo- 
meter temperature is probably real. But, prior to lunation 19 at Hadley Rille 
there is a gradual upward drift of the temperature difference between thermo- 
couple 1 and the uppermost platinum resistance thermometer and an abrupt 
decrease occurs between lunation 18 and 19. It is impossible to tell whether 
these changes are due to changes in the reference thermometer calibration, 
or the changes in EMF of the thermocouple 1 circuit. If the reference thermo- 
meter readings are corrected assuming the thermocouple drift up to lunation 
18 is entirely to reference thermometer calibration drift, then the resulting 
mean electronics box temperatures show a decrease with '.me and an abrupt 
change at lunation 19 which is unlikely. Consequently, we think that the 
reference thermometer temperatures are accurate as they are shown in 
Figure 6C, and that the characteristics of thermocouple 1 have changed. 

The presunrise temperatures of the reference thermometers are 
thermostatically controlled, and thus, not determined by the radiative balance 
of the electronics box with its environment. The differences in presunrise 
temperature between thermocouple 1 and the uppermost platinum resistance 
thermometers are stable at both experiments, which lends more support to 


I 




our assertion that the reference thermometers are stable to within a few tenths 


of a degree over three years. Although the possibility of instrument insta- 
bility cannot be ruled out entirely, for all of the reasons given above it appears 
that the upward drift of the reference thermometer temperature at noon is a 
real rise in electronics box temperature. 

This slow increase in temperature with time has been observed by 
thermometers on other ALSEP experiments. However, the results from the 
heat flow experiment are probably the most accurate. The most obvious 
cause of this temperature rise is a gradual increase in the absorptivity- 
emissivity ratio of the surfaces of the electronics box with time. This could 
result either from a degradation of the painted surfaces of the electronics 
housing or from a slow accumulation of fine lunar dust on the outer surfaces 
of the housing. We certainly can exclude the possibility that it represents 
a general rise in lunar surface temperature. 



ConcluBionB 


Equation (1), used in conjunction with a finite difference model 
the lunar regolith and the insolation function, can be used to explain the 
long term temperature variations of the exposed thermocouples at Apollo 
sites 15 and 17. These thermocouples experience strong annual variations 
i 1 temperature with amplitude modulations due to the precession of the 
;.aoon. The amplitude and phase of these modulations depend on cable 
orientation. The data presented here indicate no detectable secular drift 
in lunar surface temperature over the three year observation period. 

In situ comparison of temperatures of three different types 
of thermometers indicates that the Heat Flow Experiment thermometers 
and detection circuits are stable to within a few tenths of a degree over the 
periods of observation, except for the drift in junction 1 , probe 1, at 
Apcllo 15. The gradual upward drift in the near-noon temperatures of the 
reference thermometers indicate a steady increase in electronics box 
temperature at the two sites. A likely explanation of this increase is an 
increase in ha ratio of optical absorptivity to infrared emissivity of the 
surface* 0. the electronics housing, due either to a slow accumulation of 
duot On the housing surfaces ora gradual degradation in the reflective 
properties of the paint on the housing. 
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Figure 3 - Inset; three thermocouple temperature variations for a 
typical lunation; A-maximum and presunrise temperatures for a 
thermocouple above the surface at Taurus Littrow; B-similar data 
for two thermocouples at Hadley Rille; thermocouple 4 is located 
inside the borestem projecting above the surface. 
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Figure 4 - A; the variation in insolation period length at Hadley Rille as 
predicted by the insolation function derived in the Appendix (£q. 11) and 
compared to the observed length based on abrupt changes in thermocouple 
temperatures - B;theoretical variation in insolation intensity maxima at 
Apollo 15 and 17 sites. 
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Figure 5 - A; maximum thermocouple temperatures compared with lunar 
surface temperatures predicted by the insulation function and a thermal 
model of the lunar regolith - B; theoretical and observed presunrise 
temperatures. 
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Figure 6 - A; maximum and preaunrise reference thermometer tem- 
peratures at Taurus Littrow - B; temperature difference between 
thermocouple 1 and the uppermost platinum resistance thermometer 
in probe 1 at Apollo 17 near noon (maximum) and just before sunrise 
C and D; similar observations at Hadley Rille. 
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APPENDIX 

LUNAR INSOLATION FUNCTION 

I. Motions Considered 

A simple expression for the solar flux at any point on the surface 
of a smooth, spherical moon has been derived, based or tout idealized 
motions : 

1) Annual revolution of the earth around the sun, assumed tc be 
Keplerian and counterclockwise; 

2) Synodic revolution of the moon around the earth, assumed to 
be uniform circular and counterclockwise in the ecliptic plane; 

31 Uniform rotation of the moon upon its spin axis, counterclock' 
wise with a Draconitic period ( the interval between two successive tran- 
sitions of the moon at its ascending node); and, 

4) Uniform precession of the spin axis about a normal to the 
ecliptic plane at a constant angle of 1. 5*, clockwise with a period of 
18. 6 years. 

II. Derivation 

The flux s on a unit surface area perpendicukvr to u^ at A” (see 
{‘igure Al) is equal to the flux on a unit area normal to the sun line 
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u** which is considered to be parallel to u , times the cosine of the angle 
s s 


oc between u*' and ul* 
s A 


s = s cos oc • 
n 


s depends only on the distance of the moon from the sun. Cos oc is equal 


n 


to the pro; ction of u into u . 

A. S 


A. Determination of s . 

n 




( 1 ) 


where Q is the source strength of the sun and is the sun-moon distance 
(see Figure A2). By the law of cosines 




=Rg (1+ - 2 ) cos ♦) 

El E 


r 25 

Since — ® • 0027, the square of this ratio will be ignored, and 

R 93 


(1 - .0054 cos ) 


( 2 ) 


R can be written as 
£ 


R = a(l - e cos { )* where 
£ 


( 3 ) 


*see for example, W. M. Smarts Textbook on Spherical Astronomy 
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a is the length of the semi-major axis of the earth's orbit, e is the eccentricity 

of the earth's orbit, and { is the earth's eccentric anomaly, which can be 

expanded in terms of the eccentricity and the mean anomalyc^as 

3 Z 

{ +(e-e /8) sint^ +1/2 e sin ZcM>\ 

(for t, in mean solar d&ys.ciC - 2 (t-t )/T, and T = 365. 25 days. The 

P 

subscript p refers to perihelion). Since e is about .017 , terms involving 
higher powers of e than the first will be ignored. So 

+ e sin Mi> (4) 

combining (3) and (4) gives 

R a(l - e cos(*^+ e sincMi)) 

= a(l - e(co8<44ffcos (e sin M) - sin sin (e sin tM’)). 


Squaring both sides, noting that cos (e sinUL) ^ cos e = 0. 99986 and 
sin (e sincM^) ^ e ain cMs, and dropping the term in e^ leaves 




1 - 2e cos 


M). 


(5) 


Combining (2) and (5) gives 

R ^ a ( 1 -0, 034 cos 1 - 0. 0054 CDS ^ ) 


M 


and from (1) 


s 

n 



1 

(1-0. 034 cos4,i4)(l - 0. 0054 cos > ) 


Since 0. 034 and 0. 054 are both small compared to I, and the 
product is less than e^, s^ can be sufficiently approximated as 
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8 "= ) (1 + 0. 034 cos4/W?+ 0, 054 cos ^ ) 

n a 

using the expansion 1/1 - x= 1 4-x - x^/2 + 

2 

If we define S as Q/a , which is approximately the mean flux 
over a year (the solar constant), then 

s S(1 + 0.034 cosUt+ 0.0054 cos (6) 

n 

B. Determination of Cos oc 

As the moon rotates on its spin axis, the point A" describes 
the circle shown in Figure 1. Figure A3 shows this circle projected 
onto the ecliptic plane along 0*0**Z*. The semi-minor axis of the 
ellipse lies along O O’, and the semi-major axis along O’X’, The 
length Df the semi-major axis is sin X which is the length of the vector ^ 
The length of the semi-minor axis is sin X cos 0. Considering the sun 
to be infinitely distant, 

cos oc = (0 + ^ ) (7) 

where P denotes the length of and 0 and fi are as showii. 

Expanding (7) yields . 

cos oc = /°(cos 0 cos 0 - sin 0 sin fi ) (8) 

From Figure A3, cos 0 - x! **cos ! p - sin X cos w//^ (9) 


and sin 0 - y//^ = (Ay + yM//^ = (cos X sin 9 + sin X cos 0 sin w )//^ (10) 
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) 


By (8). (9). and (10) 

cos oc = cos ^sinXcosu- sin^cosx sinG - sin^sinX cos 9 sinu 
= sin X (cos^cos <d - cosGsin^sinw - cot X sinGsin^) 

^ sin\ {cos{<p +v») - cotX sinG sin0) , (11) 

since cos 9 (= 0.99966) is nearly 1. 

C. The Angle <f> 

0 depends on the earth's true anomaly, v, the correction to 
this anomaly for the moon's revolution about the earth, Av (see Figure A2), 
and the precession > . v can be expanded in terms ofc/i^and e as 

13 5 Z) 

V = o(4+ (2e - “ e ) sint>/^+ "re sin ZtMp + 

4 4 

2e sincJ^ (12) 

to the same approximation used for E in equation (4). 


By the law of sines (in Figure 2) 
r R. 




- sinAv sin ^ 


or sinAv = - r sin ^ /R 


M 


Substituting for R from (2) and using sinAv ^ Av, 
M 


Av - - 0. 0027 sin 4 


( 13 ) 


to the level of approximation being used. 



Letting the subscript 0 denote initial values 

+ O' 0027 7^) (14) 

using (13), with v given by (12) and where 0 and 7 are as described in I. 

Ill, Time Origin and Ir' ial Conditions 

Since the eccentric anomaly ( and the true anomaly v are 

computed from perihelion, it is convenient to use the occurrence of 

perihelion on January 2. 50, 1973 as a time origin. Thus, v = 0. 

o 

The value of can then be calculated by noting that ^ = 0 when the sub- 
solar point crosses the moon's equator going north, which occurred in 
1973 on January 8. 65, and that ^ = 0 a.t new moon which occurred on 
January 4. 67. Vcan be arbitrarily set to = 0 at perihelion 1973. 
Using values from the American Ephemeris and Nautical Almanac (14) 
becomes 

- 0.112 - v(t - t ) + 0. 0027 sin( •/» (t - t ) - 0. 46) - 7 (t - t ) 
o o o 

measured in radians with t in mean solar days. 

The initial value for oi is the difference in selenographic longi- 
tude between the point A and the sun at t = t . From the Almanac 

o 

<w = 4. 20 - col 
. o A 

where col denotes the selenographic colongitude of A. 
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IV. Comparison with Apollo Data 

The roots of the expression inside the parentheses of equation (11) 
correspond to sunrise and sunset times on a smooth, spherical moon. 
These times can be compared with sunrise and sunset times determined 
from the exposed thermocouples at the Apollo 15 and 17 sites. They 
agree within the limits of accuracy of the topographic corrections. The 
intervals between consecutive sunrises as predicted by equation (11) dif- 
fer in a random way (no systematic error) from the intervals determined 
from the thermocouple data over the times the experiments have been in 
operation with a standard deviation of 0. 028 days or about 40 minutes. 

The thermocouple data are sampled about every hour, which would give 
a standard deviation of about 0. 017 days, or about 24 minutes, assuming 
that the errors in estimation of these sunrise times are uniformly dis- 
tributed over an hour interval. The other 16 minutes are assumed to be 
due to the effects of topography and possibly to inaccuracies in the 
determination of the initial angles. 



X 


Figure A / 


The moon's spin axis (00"Z") is inclined at a constant angle 0 with rerpect 
to OZ which is normal to the ecliptic plane. 0'0"Z' is parallel to OZ, OS 
(u^), OX, O'X', OO'Y(Y'), and OA ( all lie in the ecliptic plane. The 
angles 0 and are measured in this plane. The lines OX, O'X'and 0"X" 
are considered to be mutually parallel: ~Ug and Ug" are considered coplanar 
and parallel; therefore, oc is the angle measured from vT'g to u^" alsa. (The 
sun is assumed to be infinitely distant for these approximations. ) \ is the 
colatitude of the point at A" measured from the North Pole. The u's are 
unit vectors. 




